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Abstract 

Statistical systems, in which spontaneous symmetry breaking can be accompanied 
by spontaneous local symmetry restoration, are considered. A general approach to de- 
scribing such systems is formulated, based on the notion of weighted Hilbert spaces and 
configuration averaging. The approach is illustrated by the example of a ferroelectric 
with mesoscopic fluctuations of paraelectric phase. The influence of the local symme- 
try restoration on the system characteristics, such as sound velocity and Debye- Waller 
factor, is discussed. 
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1 Introduction 



It is generally accepted that phase transitions are related to symmetry changes, so that 
the transformation from a disordered to an ordered phase is accompanied by spontaneous 
symmetry breaking [TH3]. There also exist more complicated cases, when the symmetry, 
being broken in the major part of the sample, at the same time, is restored in other parts 
of the same system, which is called symmetry restoration Reciprocally, in the generally 
disordered phase, there can appear local regions with broken symmetry. In dynamical mod- 
eling, the arising regions of symmetry that differs from the symmetry of the surrounding 
matrix are related to the generation of solitonic clusters and local coherent structures [5HIH] ■ 
It turns out that the systems with such local structures are more stable, as compared to 
homogeneous systems. Generally, the systems that become more stable by spontaneously 
changing their properties pertain to the class of self-optimizing systems [20J. There are nu- 
merous examples of condensed matter, where the sample is not homogeneous, but consists of 
regions of different symmetry. One tells that such systems display mesoscopic phase separa- 
tion, and they are termed heterophase. This, for instance, concerns many high-temperature 
superconductors, in which superconducting regions coexist with normal regions [2TH25] , and 
some low-temperature superconductors [26l[27]. Such a coexistence of superconducting and 
normal phases can occur even in atomic nuclei [28]. In many magnetic materials, magneti- 
cally ordered phase includes paramagnetic clusters [29H88] . Around the points of structural 
phase transitions, there exist regions, where phases with different symmetry coexist [3H - HU] . 
A number of ferroelectric materials displays the coexistence of ferroelectric and paraelectric 
phases [HH5l], which also concerns such novel materials as relaxor ferroelectrics [551458] . 
Much more examples of heterophase matter can be found in the review articles [qTJ|I6T)] . 

The basic difficulty in the description of heterophase systems is the necessity of dealing 
with two or more different phases, possessing principally different symmetry properties, but 
coexisting inside the same volume. The problem is aggravated by the fact that the location 
of the germs of different phases in space is chaotic. The situation is drastically different 
from the case of a sample consisting of several domains with well defined spatial locations 
and structure [SI]. A typical heterophase system is a mixture of regions, randomly located 
in space and having various and often very irregular shapes. Moreover, in many cases, the 
heterophase regions are not static, but can move in space, vary their shapes, and even appear 
and disappear. 

To describe such a complicated matter, it has been necessary to develop an approach 
allowing for the treatment of these inhomogeneous and nonequilibrium systems. More pre- 
cisely, such systems are to be locally equilibrium, since the notion of phase requires the 
existence of at least local equilibrium. A general approach for treating heterophase systems 
has been advanced [62H66] and summarized in reviews [59|. [60] . 

One of the main problems in treating the systems with coexisting regions of different 
symmetry is how to separate the states, corresponding to different symmetry properties, 
in the space of microscopic quantum states. The standard situation is when the system 
as a whole is characterized by a given Hilbert space, with a prescribed symmetry. Then 
how would it be possible to describe subsystems of different symmetry in the frame of the 
same Hilbert space? This can be done, for instance, by imposing different conditions on 
the equations characterizing different phases [H2JE7] or by invoking the method of restricted 
trace [SS]. Probably, the most powerful, convenient, and rigorous method of separating 
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different phases is the approach based on the notion of weighted Hilbert spaces, advanced in 
Refs. [59(160] . In the latter publications, however, this method was not formulated in the most 
mathematically general way. The aim of the present paper is to further develop the method 
of weighted Hilbert spaces, formulating it in the most general form. As an illustration, 
the method will be applied to deriving a general model of a heterophase ferroelectric. The 
choice of ferroelectrics for illustration is caused by several reasons. First, there are numerous 
examples of the heterophase systems of this type [H1460] . hence their correct description is of 
great importance. Second, the general model related to these materials is generic for many 
other substances, hence it serves as a good example for extending the approach to other 
heterophase systems. 



2 Single-Phase Systems 

Before developing an approach for treating multiphase systems, it is necessary to briefly 
recall the main points of describing the standard case of a single-phase system, which will be 
used in what follows. For generality, quantum systems are considered throughout the paper. 
Let us have a set {tp n } of states forming a basis for the closed linear envelope 

£ = Span n {(^ n } . (1) 

Here the index n implies a multi-index that can pertain to either discrete or continuous set. 
The discrete set can be infinite. Let the scalar product < f\h > be defined for each pair 
f,h £ £. The norm, generated by the scalar product, is 

11/11 = >/<7T7> Cfe£)- (2) 

Completing the linear envelope ([[]) by the norm yields the complete normed space, that is, 
the Hilbert space 

n = {£, ||/||} . (3) 

The so-defined Hilbert space can be separable or not, depending on the physics of a concrete 
problem and, respectively, on the nature of the multi-index n enumerating the basis {f n }- 
The existence of a Hilbert space, associated with the considered physical system, is the 
necessary prerequisite for characterizing the system. 

The basis {(p n } can be taken to be ortho normalized, such that 

( V?m I ) = S mn . (4) 

And let it be complete. Then for each f,hES, one can write their expansions over the 
basis, 

f =/; fnVn , h = ^ h n¥n , (5) 

n n 

with the expansion coefficients 

fn = ( | / ) , K= (fn\h) . 

Therefore, the scalar product can be represented as 

(f\h)=J2 fnK- (6) 
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Suppose the algebra A of local observables, represented by Hermitian operators A e A, 
be given on %. The matrix elements of A over the basis {p n } are 

Arm = ((p m \ A \ (p n ) . (7) 

The system statistics are characterized by a statistical operator p that is a trace-one operator, 
acting on %. The operator averages are defined as 

(A) = Tt h pA = ( | f>A | <p n ) . (8) 

n 

The set of all available averages is termed the statistical state: 

{A) = {(A)} . (9) 

Characterizing thermodynamics phases, a special role is played by the order operator 
fj £ A, whose average defines the system order parameter 

v=(v). (10) 

In general, this can be a scalar, vector, or matrix quantity, so that one talks of an order 
parameter just for short. The main point is that this order parameter is specific for each 
thermodynamic phase, such that different phases, possessing different symmetries, enjoy 
different order parameters. 



3 Symmetry Breaking 

In the language of symmetries, the existence of phase transitions assumes the following. If 
there is a statistical state that is invariant with respect to a symmetry group, then it can be 
decomposed into a sum of several terms describing different pure phases [U[69l[70] . To select 
a state with a particular type of symmetry, one employs one of the variants of the quasiaver- 
aging techniques [UEZ1EE] • When the state with a broken symmetry is thermodynamically 
more stable than the invariant state, one calls this the spontaneous breaking of symmetry. 

The methods of quasiaverages allow one to select a particular state with the desired 
symmetry only in the case of thermodynamically equilibrium systems, when all the system 
is characterized by one and the same symmetry. But our aim is to describe the situation, 
when inside the considered system there appear regions with different types of symmetry. 
How could we proceed in such a principally different case? For this purpose, it is necessary 
to develop a more general method of symmetry breaking, which could be used not only 
for equilibrium systems, but also for quasiequilibrium, metastable, or even for arbitrary 
nonequilibrium systems. Such a general method is developed below. 

Let us consider a statistical system that, under different thermodynamic conditions, could 
be in different thermodynamic phases, enumerated by the index v = 1,2,.... We keep in 
mind that the space of microscopic states, related to the considered system, is the Hilbert 
space (j3J), with a basis {<p n }, as described in Sec. 2. Let us put into correspondence to a 
vector p n a probability 

Pn=P U M- (11) 
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The set {p 1 ^} of these probabilities forms a probability measure with the standard normal- 
ization property 

E Pn = 1 (0 < K < 1) ■ (12) 



Now, let us introduce the weighting operator 

n 

Using this, we define the weighted scalar product 

(f\h) v = (f\W v \h) . 
With the weighting operator ([TBI , we have 

(f\h)u = '£ l ft(f\<Pn)(<Pn\h). 
n 

Under expansions the latter reads as 

n 

The scalar product (TL5I) generates the norm 

ll/ll, = >/</!/>, (/eS), 

which, in view of form ffl6|) . can be written as 



(13) 



(14) 



(15) 



(16) 



(17) 



'£ Pn W ■ (18) 

n 

The closed linear envelope ([TJ, equipped with norm ( 1171) . is the weighted Hilbert space 

n v = {£, 11/11,} . (19) 



The representation ^ of the algebra of local observables A, acting on the weighted 
Hilbert space ffl9|) , consists of the operators A„ defined through their matrix elements 



Taking into account the weighting operator ( fl~3|) gives 
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+ 



( <^ m | a v | v? n ) = - ( p v m + pj; ) ( ip m I a i <p n ) . 



The operator averages are 

(A v ) = Tl Uu ( pA ) = J] ( V» I (M)^ I V?n ) 



(20) 



(21) 



(22) 



With the matrix elements (j2ip . we get 

(A„) =J2Pn(V>n\pA\ip n ) . (23) 
n 

Similarly, the average of the order operator fj u G A v becomes 
Respectively, the order parameter is 

^ = ( Vu ) ■ (25) 

The probability measure {p^} is to be defined so that to guarantee the order parameter, 
specifying the chosen thermodynamics phase. This means that the probabilities ( ITTj) should 
select the states typical of the considered phase [72J[73], which does not need to be equilib- 
rium. By this construction, it is clear that the scheme, based on the weighted Hilbert spaces, 
includes as a particular case the selection of phases by means of the quasiaveraging method, 
since the latter also chooses the states typical of the desired phase, but provided this phase 
corresponds to a stable equilibrium system. 

4 Multiphase Systems 

First of all, let us stress the difference of the case we try to describe, as compared to the 
Gibbs phase mixture [73]. In the latter case, the system is spatially separated into several 
macroscopic regions filled by different thermodynamics phases. But in the case we are 
interested in, the system is a heterophase mixture, comprising the mesoscopic germs of 
several thermodynamic phases, which are randomly intermixed and coexist in a region of 
thermodynamic parameters in the same volume [591160] . The space of states for such a 
heterophase system is the fiber space 

7 = ^U U . (26) 
The basis of the latter, {</? n }, is made of the tensor products 

£n = (g) <fin„ , (27) 
v 

in which the notation for the multi-index 

n = {m, n 2 , n 3 , . . .} 
is used. Any state / G J can be decomposed over this basis: 

7= Yl f n ^ n ( /n = ( | 7) ) • 
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Then the scalar product of any pair /, h 6 J 7 is given by 

{J\h)=Y,^K- (2* 

n 

The scalar product generates the norm 



(/!/>■ (29) 
The fiber space (j26j) can be equivalently represented as the closed linear envelope 

^={Span n {^ n }, 1 17| | } , (30) 



equipped with norm ( 12911 . 

The operators of observables A, acting on J 7 , are defined as the direct sums 



A = K • (31) 

V 

The related matrix elements are given by the expressions 

( £ m | A | (p n ) = ( <p mv \ A v | y nv ) Yl ^m M n M ■ (32) 

The operator averages are defined as 

(A) = Tr^ pA = J2 (<P* I M I ) • (33) 

n 

This, in view of form ( l27l) . yields 

&) = '$2'52(<Pn\{pA) v \<p n ). (34) 



Comparing this with ( l22j) . we obtain 

(I) = ^ (I,) . (35) 

Z' 

The set {< A >} of all observable averages is the statistical state of the heterophase system. 

5 Phase Configurations 

When the system is inhomogeneous, being composed of many mesoscopic regions of different 
phases, we need, first of all, to describe the spatial distribution of these regions inside the 
system. For this purpose, the space V, occupied by the system, can be decomposed into the 
subregions, whose set {V„} forms an orthogonal covering: 
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such that 

v„ n v " = ' ( 3? ) 

where 

V = / dr , V u = f dr. (38) 



The set of the regions, occupied by a //-phase is described by the manifold characteristic 
function [75], or manifold indicator 

«-« s {i: r^vl • < 39 > 

These functions satisfy the properties 

]T £„(r) = 1 (r e V) (40) 

and 

r &,(r) dr = K • (41) 



The collection of all manifold indicators defines the phase configuration 

£ = {&,(r) : !/ = 1, 2, . . . ; r e V} . (42) 

In turn, the space, occupied by a //-phase, can be decomposed into subregions V„, forming 
an orthogonal subcovering {V u }, such that 



V„ = Qv*, (43) 



i=l 

and 



V„i f| V ^ = MtfV* • (44) 
Each subregion Y ui is described by its manifold indicator 



in which a^i is a fixed vector pertaining to V„j. Then the manifold indicator fl39l) writes as 

u*) = Y, ■ ( 46 ) 



i=l 



Under a given configuration, the relative volume, occupied by a //-phase, is characterized 
by its geometrical fraction 

x ^vl Ur)dr= v- < 47) 

As is clear, the latter satisfy the normalization condition 

J2 X » = 1 (0 <x v <\). (48) 



This shows that the set {x u } of all admissible geometrical fractions, with the given normal- 
ization condition, can be regarded as a probability measure. 

Since the the phase regions are randomly distributed inside the considered system, the 
configuration f l4"2"j) is to be treated as a random variable. The related differential measure 
can be defined as 

^ = 'fc^-'ln^nfiT' (49) 

\ U / V V 1=1 

All possible configurations constitute a topological space 

X = {£, V$ . (50) 

The above constructions describe the situation, when the system is separated into several 
regions filled by different phases, and these regions are randomly distributed in space. 



6 Configuration Averaging 

An inhomogeneous heterophase system is not necessarily in complete equilibrium. But it 
must be at least quasi-equilibrium in order that it would be admissible to talk about the 
germs of phases. As is known, the notion of phases is not strictly defined for finite systems. 
The mathematically rigorous definition of phases assumes the introduction of thermodynamic 
limit [2l[3],[72l[7U[76]. However, in practice, it is possible to speak about germs of phases 
already when each of such germs consists of a large number of particles iV ^> 1. Monte 
Carlo simulations show that thermodynamics phases can be well defined already for 10 — 100 
particles in a finite cluster [77] . Respectively, though the symmetry, related to a phase, may 
be not strictly defined for a finite cluster, but it is possible to talk about an asymptotic 
symmetry that is approximately defined for a large number of particles N 3> 1, keeping in 
mind that the symmetry becomes exact in the thermodynamic limit. 

A quasi-equilibrium system is described by a quasi- Hamiltonian [59] depending on a given 
phase configuration and having the operator structure as in (13T|) : 

QiA) = @QuiAu). (si) 

V 

The partition function 

Z = Tr^J exp{-Q(£)} P£ (52) 

includes the quantum averaging over the given quantum variables and the averaging over 
phase configurations. This function defines the quasi-equilibrium thermodynamic potential 

y = -l\nZ. (53) 

Considering asymptotically large systems, we, as usual, keep in mind the thermodynamic 
limit 

N 

N ->■ oo , V — > oo , — ->■ const . (54) 
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Similarly to (131]) . the operators of observables, under the given phase configuration, have 
the form 

i(0 = 04,(60- (55) 

And for what follows, it is convenient to introduce the notation 

A v {x v ) = lim A V {Q . (56) 

The statistical operator of a multi-phase system, with a fixed phase configuration, is 

P(0 = 4 exp{-g(0) . (57) 



z 

The observable quantities, related to the operators (155]) . are given by the averages 

(l)^Tr^ | p(£)A(Z)V£. (58) 

The following theorem is valid [59 ] 178 ] 179] . 
Theorem 1. 

If Qu{^ u ) can be represented as an expansion in powers of £„, then the thermodynamic 
potential (153]) . in the thermodynamic limit ( 15%]) . asymptotically equals 

y = abs miny(w) , (59) 

w 

with the set 

w = {W!, w 2 , ■ ■ .} (60) 
forming the probability measure enjoying the standard properties 

J2w u = l (0 < w„ < 1) , (61) 



and where 



is the sum of the terms 



in which 



y{w) = ^ y u (w u ) (62) 



Vv{wv) = - InTr Z u , (63) 
Z v = Tr Uv exp | -Q„(w v )j . (64) 



The quantities w v are the phase geometric probabilities, showing the corresponding weights 
of the coexisting thermodynamic phases. 

Corollary 1 
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From this theorem, it follows that the thermodynamic potential fl59|) can also be repre- 
sented in the form ^ 

y = - — In Tr Z , (65) 

where 

z=j[z v . 

V 

The above theorem defines the thermodynamic potential of a heterophase system. The 
observable quantities, corresponding to the averages of the operators from the algebra of 
local observables, are given by the following theorem [59 t l78 | [79]. 

Theorem 2. 

Assume that Q u (£ u ) and A u (£ u ) can be expanded in powers of £„, then the averages f )58|) . 
for asymptotically large N, take the form of the sum 

(A) = £ A) , (66) 

with the terms 

(A u ) = Tr w „ p v A v , (67) 
in which ^ 

^ = ~z~ exp ' ( 68 ) 

and where the notation A v = is used. 

Corollary 2. 

Equivalently, the averages ( 166]) can be represented in another way by introducing the 
operators 

A = @A V , (69) 

V 

for which 

(A) = Tr f pA , (70) 

where 



V 



The second theorem defines the method of calculating the averages for a heterophase 
system. 

7 Effective Hamiltonians 

The quasi-Hamiltonian Qv{£v) is connected with the local Hamiltonian H u (r,^ u ) through 
the relation 



Qu{£u)= / Mr^ v )H v {v,i v ) dv , (71) 
11 



in which the inverse local temperatures (3 v (r, £ u ) play the role of the Lagrange multipliers. 

If the considered system is such that the phase fluctuations can appear randomly in any 
part of the system, the latter is called phase-uniform on average [59] . For such a system, the 
configuration averaging of the inverse temperatures gives 

Pu = J P v (t,Z„)VZ. (72) 

The mixture of different phases is in thermal equilibrium, when the temperatures of these 
phases coincide: 

& = V = \ W ■ (73) 
Then the quasi-Hamiltonian Q v (w v ), entering the partition function fl64|) . becomes 

Q v (w v ) = V , (74) 

with the renormalized Hamiltonian 

Hu = J H u (r,w u ) dr . (75) 
The partial statistical operators ( 168)) acquire the form 

Pv = y exp (-pfiSj , (76) 

with the partition functions 

Z„ = Tr n „ exp . (77) 

It is convenient to introduce the effective Hamiltonian 

H = Q)H U , (78) 



using which, the thermodynamic potential ( !65j) can be represented as 

y = - 1 In Tr e"^ . (79) 

The terms of sum (1781) are called the phase-replica Hamiltonians, since they have a similar 
mathematical structure, but are associated with different phases. 

In order to connect the thermodynamic potential ( 1791 with the free energy, let us recall 
that the latter is defined as 

F = —Tin Tr e~ pH . (80) 
Therefore, the potential y is directly connected to the free energy by means of the relations 

F = NTy = Y J F v 

and 

F, = -Tln Tr Z v = NTy u . 
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To specify the consideration, let us take the local Hamiltonians in the usual form 



Hv{r,t») = &(r)^(r) (-^- + U) ^(r) + 

+ \f Ur)Ur')ri(rM(r'Mr ~ OW^M dr . (81) 

Here U — U(r) is an external potential and ?/v( r ) are the field operators of the particles 
forming the system. 

Then the phase-replica Hamiltonians are 

Hu = w v J 4>l{r) (-¥. + u\ ^(r) dr + 

+ Yj ^M^Mr - v')^{v')^{v) dvdv' . (82) 
Introducing the notation for the kinetic-energy operator 

K v = [ $(r) f-^l + u] ^,(r) dr (83) 



2m 

and for the operator 

$ v = J ^(r)^( r ')$( r _ rOVv(r')Vv(r) drdr' , (84) 
related to the potential energy part of the Hamiltonian ( )82|) . we have for the latter 

2 

H v = w u K u + ^r®v (85) 

In this way, we have derived all basic equations for treating heterophase systems. The 
derivation has been based on the following three major points making it possible to separate 
different thermodynamic phases. First, to distinguish the phases in the space of microscopic 
states, the notion of weighted Hilbert spaces is introduced. Second, to separate the phases in 
real space, the manifold indicators were employed. And, finally, the procedure of averaging 
over phase configurations is accomplished, leading to the set of equations for equilibrium 
on average phase replicas. The idea of the averaging procedure reminds the method of 
averaging [80] and the scale separation approach [81TI83] . used for nonlinear equations in 
dynamical theory. The main difference from the latter is that here we have averaged out 
slow heterophase fluctuations, slow with respect to the fast microscopic motion of particles, 
while in dynamical theory one usually averages out fast fluctuations, leaving at the end the 
slow motion of guiding centers. 



8 Stability Conditions 

The developed theory should be complimented by an important addition discussing the sta- 
bility of heterophase systems. Thermodynamic stability is characterized by the minimization 
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of a thermodynamic potential. For concreteness, let us consider the case of a heterophase 
system, where two phases coexist, so that v — 1, 2. And let us denote 

w\ = w , w 2 = I — w . (86) 

Minimizing the thermodynamic potential f!79|) with respect to w implies 



dy d 2 y ^ 

<9u> <9u> 2 



The first of these equations gives 

<£>-»■ 

And the second condition leads to the inequality 

, <9 2 # . nl fdH\\ , , 

<^> >,? <(^)>- < 89 > 

Since the right-hand side in the above inequality is nonnegative, the necessary condition of 
heterophase stability is 

To specify these conditions, let us take the effective Hamiltonian (ITS]) , with the replica 
Hamiltonians ( )85l) . And let us use the notations for the averages 



K v EE , $„ EE ($„) . (91) 

Then we have _ 

1 <9i^ 

— ( ) = K x + w$i - K 2 - (1 - u;)$ 2 
TV 

and _ 

This yields the equation for the phase probability 

$ 2 + K 2 - 

w = ^ t • (92) 

$1 + $2 

The stability condition f[8"9~j) gives 

>l<(f)>< < 93 » 

and from the stability condition fl90|) . we get 

$i + $ 2 > . (94) 
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One more condition, follows from the definition of w„ as of phase probabilities, according 
to which 

< w < 1 . 

This, with the use of (1921) . results in the inequalities 

- $1 <Kx-K 2 < $ 2 • (95) 

The stability conditions must be valid in order that the considered heterophase system 
could really exist. Note that the stability condition ( 1901) is analogous to the condition of 
diffusion stability [83] . 



9 Heterophase Ferroelectrics 

To illustrate how the developed theory works, let us consider ferroelectric materials. Ferro- 
electrics are known to be a good example of matter demonstrating heterophase properties. 
In many ferroelectrics, above the transition point from the disordered into the ordered state, 
there exist small polarized clusters [37] . These fluctuational embryos of the ordered phase 
inside the disordered phase where termed by Cook [36] antiphase fluctuations, emphasizing 
that they where a particular case of heterophase fluctuations, whose existence was predicted 
by Frenkel [85]. Such embryos of one phase inside another phase can exist in a whole re- 
gion around the phase-transition point T c . At temperatures below T c , the embryos of the 
paraelectric phase inside the ferroelectric phase arise at a temperature T n called the lower 
nucleation point [SB]. And above T c , there is another temperature T*, called the upper 
nucleation point [SB], where the embryos of the ferroelectric phase appear inside the para- 
electric phase. Thus, around the phase transition temperature, there can exist a region of 
temperatures 

T < T < T* 

where the embryos of two phases coexist. The phase transition, generically, can be either 
of first or of second order, but the appearance of heterophase fluctuations, usually smears it 
into a continuous crossover [4TH33]. 

Heterophase fluctuations were observed, for instance, in such well known ferroelectrics as 
HC1, DC1, mixed crystals HCli_ x DCl x , and RbCaF3, where they were intensively studied by 
Brookeman and Rigamonti using nuclear quadrupole resonance [4~2ll4"3] and nuclear magnetic 
resonance [2] . Heterophase fluctuations in these ferroelectrics arise in a finite region around 
the phase transition point. The appearance of these fluctuations occurs even without external 
defects, although the presence of defects intensifies their nucleation [JT]. 

Such heterophase fluctuations have also been observed in many other ferroelectrics, e.g., 
in C 4 4 H 2 [15, in KH 2 As 4 @6], in Rb^ND^i^DaPC^ [S], in Na x Bi!_ x Ti0 3 [49], and 
others. They also arise in such novel materials as relaxor ferroelectrics, for example, in 
(PbZn^Nbi^Og) [53], in Ba 2 NdTi 2 Nb 3 0i 5 and Ba 2 La . 5 Ndo.5Ti 2 Nb 3 Oi 5 [SgJES], and, with a 
high probability, in many other relaxors, such as PbMgi/3Nb 2 /303, (Ba x Pbi_ z )(Zni/3Nb 2 /3)03, 
(Sr 2: Pb 1 -. :c )(Zn 1/ 3Nb 2/ 3)03, (Ba x Pbi_ x )(Ybo. 5 Nbo.5)0 3 , Pb 1 _ z Ba :c (Yb 1/2 Ta 1/2 )03, (Fe 1/2 Nb 1/2 ) 
Pb(F ei / 2 Ta 1/2 )0 3 , Pb(Yb 1/2 Ta 1/2 )0 3 -Pb(Fe 1/2 Ta 1/2 )0 3 , and Pb(Mg 1/3 Nb 2/ 3)0 3 -PbTi03, [58]. 

A simple model of a heterophase ferroelectric has been considered in Refs. [5T H IS7 ] ISS]. 
In this model, phonon degrees of freedom were not taken into account. However, the latter 
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are important because of two reasons. First, the appearance of heterophase fluctuations is 
frequently accompanied by heterostructural fluctuations [89,90j, which are connected with 
phonon excitations. Second, the occurrence of the paraelectric-ferroelectric phase transition 
is intimately related with phonon characteristics that can be directly measured. In the 
present paper, we derive a generalized model of a heterophase ferroelectric, taking into 
account the phonon degrees of freedom. This makes it possible to find out the influence of 
the heterophase fluctuations on such observable quantities as the Debye- Waller factor and 
sound velocity. This influence is especially pronounced in the phase-transition region. 



10 Basic Hamiltonian 

Let us consider a ferroelectric, in which there can arise the embryos of the competing phase. 
So that the sample can house two coexisting phases, ferroelectric and paraelectric, being 
randomly intermixed with each other. The ferroelectric phase will be indexed by v = 1 and 
the paraelectric phase will be labeled by the index v — 2. For concreteness, we shall study 
the model of a ferroelectric, in which the order is characterized by the pseudospin operator 
describing the shift of a charged particle into one of the wells of a double-well potential at 
the cite j of the crystalline lattice [9T],[92]. Then the ordered ferroelectric phase corresponds 
to the nonzero order parameter 

(S|i> ^ . (96) 

On the contrary, the disordered paraelectric phase is characterized by the zero order param- 
eter 

(Sh) = • (97) 

Starting with the ferroelectric Hamiltonian, having the mathematical structure character- 
ized by the pseudospin variables j9U[92] , we follow the general scheme described above, and 
after averaging over the random phase configurations, we come to the effective Hamiltonian 
(J7HD with the replica Hamiltonians 

2 2 

3 i^j 3 

+ wiJ2 B^)SIS* U - wlJ2 K^SIS^ . (98) 

Here, the first term represents kinetic energy, A(r),B(r), and J(r) are particle interactions, 
Q is the tunneling frequency, and the abbreviated notation 

Tij =r t - Tj . (99) 

is employed. 

The phonon variables can be introduced in the standard way by defining the deviation 
uj from the lattice site with the lattice vector aj as 

Tj = SLj + Uj . (100) 

The lattice vectors are assumed to form an equilibrium lattice, being defined as the averages 

a, = ( Tj ) . (101) 
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Hence the average deviation, by definition, is zero, 

(uj> = . (102) 

The interactions, as usual, are supposed to be symmetric with respect to the spatial 
coordinate inversion: 

A(-r l3 ) = A(r t] ) , B(-v l3 ) = B(r y ) , I{-r i3 ) = J(r y ) . (103) 

In what follows, we shall also use the short-hand notation for the vector differences 

ELy = SLi - 3Lj . (104) 

and 

Uij = ui - uj . (105) 
The interactions are expanded in powers of the deviations as 

A{v i3 ) = Aij + A >% + \ E <<4 > ( 106 ) 

where 

A l3 = A&j) , A% ee , A% ee . (107) 



% — % 



The same expansions are made for B(r) and I(r). Owing to the symmetry properties (11031) . 
we have 

A.. — A .. A a — _ A a A<*P _ A^ a _ A^ a _ A a @ (1C\R\ 

A %3 — A J% , A^ — Ji^ , A^ — A i3 - — A- { — Aj { . l-LUOj 

In the Hamiltonian (1981) . the double summation over the lattice excludes the self-action 
terms with % — j. In order to simplify the notation, we can sum over all lattice sites, setting 
the diagonal elements 

An = B u = lu = 0. (109) 
The lattice is treated as ideal, because of which 

A = A ij = const ( 110 ) 

j 

does not depend on the index i. Using the ideality of the lattice, we get 

j 1 

and 

?^ = ^ = °- (112) 

Consequently, 

E^<-=°- ( 113 ) 
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The same type of expansions is accomplished for both the ferroelectric and paraelectric 
phases. Therefore, in what follows, we shall consider the resulting transformations only for 
the ferroelectric phase, keeping in mind that the same is done for the paraelectric phase. And 
to simplify the notation, we shall not write explicitly the index v — 1. Then, substituting 
the above expansions into Hamiltonian (|98| . with v — 1, and invoking the notation 

ij — * j ' 

we obtain the Hamiltonian 

A = »E ii + T NA + tEE^« - «nE*5 + 

+ y: Ui + e sx- + ^ e ) ^ 

ij \ a a/3 / 

in which w = W\. 



(114) 



(115) 



11 Pseudospin-Phonon Decoupling 

The obtained Hamiltonian is yet too complicated to be treated, and some approximation is 
required. For any approximation, we have to keep in mind that 

( u% } = <<>-< > = . (116) 

Because the phonon and pseudospin operators are of different nature, it is reasonable to 
decouple them in the second-order, with respect to the deviations, terms as 

<A 3 s l = ( u H ) s l + << (si) - ( << > < <sj- > ■ (H7) 

At the same time, the terms linear in the deviations can be left for a while, since later they 
can be dealt with by using canonical transformations. 
Using again the lattice ideality, we see that 

E«< = 2 E«<- ( 118 ) 

ij ^ 

The average < Sij > depends on the difference a^. Thence 

EtfW = a^E'«W = o. (ii9) 
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Similarly, 



■4 



E B t<< ( s f,) = 2 E B S" ( S S > <• 
Y. l t u %< w = 2 E J ? ( ^ ) << • ( 12 °) 



«? 1-3 

Let us introduce the notation for the dynamical matrix 

®f = Af + 2Bf(S* j ) - 2Jjf <S£> (121) 
and for the renormalized interactions 

By = £ y + <<«;>■ (122) 

a/3 

and 

la = h + ( 123 ) 

a/3 

Separating the non-operator energy part 

E 1 = ^NA + ^ 2 ££ ( J ? ( ^ ) - < W )) ( ) > ( 124 ) 

»j a/3 

we reduce Hamiltonian f II 1 5 It to the sum 

H\ — Ei + H ph + // ps + Hii n . (125) 
The second term here is the effective phonon Hamiltonian 

^ = »Ef| + T (126) 

j ij a/3 

The third term is the effective pseudospin Hamiltonian 

h ps = „•<>»;;• + ™ 2 £ (4-^ - 4^) ■ 0-><) 

And the last term is the linear pseudospin-phonon interaction Hamiltonian. The latter is 
obtained by invoking the properties 



E B ij U ij S ij 2 E B ij S ij U i ~ 2 X] />> '.' ,S '! 

ij ij 

£ = 2 £ W«? = - 2 E r >) s u"i ( 12 §) 



y y y 



and denoting 
which yields 



k% = irs;: - i«st 3 , (129) 



Thus, the Hamiltonian parameters are renormalized due to the interactions between pseu- 
dospins and phonons. 
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12 Dressed Phonons 

The quantization of the phonon degrees of freedom can be done in the way, similar to how it 
is done for pure crystalline phases. However, here we have to be careful, taking into account 
the presence of the factors w = Wi, characterizing the weight of the related phase. Thus, the 
eigenproblem equation for the phonon frequencies and polarization vectors takes the form 



£ ^^(ik.^ej,^^. (131) 



W 

m 



The phonon frequencies and polarization vectors can be chosen to be symmetric with respect 
to the momentum inversion, 

(d—ks LOks i ^—ks &ks • 

The polarization vectors are orthonormalized, such that 

e ks ' e ks' = $ss' i ^ e ks e ks = 'W ? 

s 

with the momentum summation over the Brillouin zone. The eigenproblem (11311) can be 
rewritten as 

' = (132) 

P 

where 

3 

The phonon quantization, in the presence of the linear pseudospin-phonon interactions, 
differs from the standard case by the necessity to involve a nonuniform canonical transfor- 
mation 

= - ^= E \/j Uks eks ( bks - b -ks) e ' k ' aj > 

% = v, + -±= £ ~1T= (** + b - k s) ^ • (133) 

The nonuniformity comes through an additional term in the expression for Uj. 
Hamiltonian ( I125P transforms to 

H 1 = E 1 + H' ph + H ps + H' lin . (134) 

The effective phonon Hamiltonian is 

H'ph = w^u ks \ b\ s b ks + ^ ) 

ks ^ ' 

with the phonon frequency given by the equation 

o w 



(135) 



^ks — — 

m 

j a/3 



£ £ ^ ■ ( 136 ) 
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The latter equation can also be written as 

^k e ks e 



W 

m 



at 



The momentum k pertains to the Brillouin zone. 

After transformation f 1 1 3 3 j) . the term, remaining from the renormalized linear pseudospin- 
phonon interaction, reads as 

4 = -2» 2 EE K Z v h ( 137 ) 

ij a 

with 

7;? = 

2N 

where the notation 

^^^exp^k-a^) (139) 

ks 

is used. The latter quantity possesses the properties 



ij P 



aB aB Ba Ba 



Combining (I137P and fl 138[) gives 



w 2 



tj fg aB 

which shows that this is an effective four-pseudospin interaction. 

In the summation over momenta in f)139p . the main contribution comes from the term 
with k = because of the fast oscillations of the exponential. Therefore, expression (11391) 
can be well approximated as 

y°f 7 a/3 = Aw Y . (141) 



ks 



Then H138[) becomes 



Owing to the property 



we have 



«/ = ^EE^4- ( 142 ) 

ij 8 



E< = -E4=» 

ij ij 



H' lm = , v« = 0. (143) 
Consequently, the term H' hn in Hamiltonian ( I134j) can be omitted. 
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This is in agreement with the following. The existence of the linear, in the deviations, 
terms in Hamiltonian (I130p . formally, leads to the fact that, according to (I133p . the av- 
erage < Uj > could be nonzero, which, however would be in contradiction with condition 
J2"|) . Therefore, the linear in deviations terms should be zero, as in Eqs. (|143|) . Also, 



the appearance of the linear terms in Hamiltonian (11301) breaks the symmetry of the initial 
Hamiltonian with respect to the inversion Uj — > —Uj. Generally, such linear terms should 
either be zero, due to symmetry properties, or are to be canceled by counterterms preserving 
the equilibrium condition (I102p . 

The same conclusion could be obtained, if, when decoupling the pseudospin and phonon 
degrees of freedom, we would employ the decoupling 

= ( «§ >Sg + u%( Sg ) - ( u% )( Sg ) , (144) 
which, in view of (I116p . yields 



Then, employing the properties 



d 

Y B ij ( S ij) = Y Bi i ( S ij > = ' 

3 1 3 

E J 5<^> = ^E^(^) = ( 146 ) 

3 1 3 

and _ 

E { S o ) 11 • E^^>= ' HIT) 

following from the ideality of the lattice, results in the Hamiltonian 

H 1 = E 1 + H ph + H ps . (148) 

The first term here is the nonoperator part HI 24)) . the second term is the effective phonon 
Hamiltonian 

Hph = w E u ks I b\ s b ks + - J . (149) 

ks ^ ' 

And the effective pseudospin Hamiltonian is 

Hp, = -wQ S J + w2 Y B v S i S 1 ~ w2 Y ^ S * S 3 ■ ( 15 °) 

The pseudospin interactions are renormalized by the existence of the phonon vibrations. 
And the phonon characteristics are renormalized due to the phonon interactions with pseu- 
dospins. In addition, all quantities are renormalized by the presence of the heterophase fluc- 
tuations. It is, therefore, possible to call the resulting effective phonons as dressed phonons. 
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13 Heterophase Fluctuations 

With the derived effective Hamiltonian ( 11481) . we can explicitly calculate all phonon charac- 
teristics. For instance, the phonon distribution is 

n ks = ( b[ s b ks ) = exp - 1 . (151) 

The deviation-deviation correlation function read as 

ks 

The mean kinetic energy becomes 

( — ) = 7^7 V^coth(^) . (153) 
x 2m 1 AN ^ \ 2T J v ' 

ks 

In order to treat the pseudospin variables, we can use the mean-field approximation for 
the products SfS^. This is possible because the critical region in ferroelectrics is known 
[91~||92] to be narrow due to the smallness of the Ginzburg number. Then, in the mean-field 
approximation for the pseudospins, we find the averages for the x-component 

(S?) = w-^r tanh (^] , (154) 

for the ^/-component 
and for the z-component 



2Hj \ 2T 



S? > = , (155) 



<s;> = w2 ( s Dj^ tanh (§) ■ ( 156 ) 

Here, the notations are used for the effective tunneling frequency 

Qj = Q- 2wB ( S? ) (157) 

and the effective field 

Hj = w yJtf + Aw 2 ! 2 ( sj ) , (158) 

in which 

B = J2 By , I = h ■ (159) 

j j 

To deal further with considering the properties of the heterophase ferroelectric, we need 
to restore the phase indices v — 1,2. The effective Hamiltonian of the heterophase system is 

H = H x H 2 . (160) 

The conditions, distinguishing the ordered and disordered phases are 

( } ^ , ( S* 2 ) = . 

All necessary equations for the ferroelectric phase are written above. The equations for the 
paraelectric phase can be obtained from the above expressions by setting there the order 
parameter < S% >= 0. The equations for the phase probabilities are given in Sec. 8. 

23 



14 Phase Transition 



The total system of equations, defining the heterophase ferroelectric, can be solved numer- 
ically. Below, we present some of the most interesting conclusions describing the influence 
of the heterophase fluctuations on the system properties. The strongest influence of these 
fluctuations occurs in the vicinity of the phase transition point, that is represented by the 
temperature 

~ (1 -b)u 

4artanh(2£) ' 1 J 

where _ _ 

~ B n 



I+B I+B 
The existence of the heterophase fluctuations can be noticed and their influence measured 
by studying, e.g., the Debye- Waller factor f DW and sound velocity s. It is convenient to 
characterize the influence of the heterophase fluctuations by comparing their values in the 
presence of the latter with the related values without these fluctuations. For example, 
comparing the sound velocity s in the heterophase system with the sound velocity So in a 
pure system without such fluctuations, it is useful to introduce the relative sound-velocity 
decrease 

Ss = . (162) 

so 

Another useful characteristic is the Debye- Waller factor /uvy that can be measured by 
x-ray scattering, coherent neutron scattering, and by Mossbauer spectroscopy. All details 
and definitions can be found in the books [93H95] . We need to compare the Debye- Waller 
factor f DW for the heterophase system and its value few f° r a pure system. Again, it is 
useful to employ the relative Debye- Waller factor decrease 

rn foW - fDW f. ao s 

OfDW = 7 ■ ( 163 ) 

JDW 

It is interesting that, as follows from numerical calculations, these quantities, the relative 
sound velocity decrease and the relative Debye- Waller factor decrease, are universal, weakly 
depending on the considered materials. For the relative sound velocity decrease at the critical 
transition point, we get 

Ss w -0.3 (T = T c ) (164) 
and for the relative Debye- Waller factor decrease at this point, 

Sfow ~ "0.3 . (165) 

This decrease of the sound velocity and of the Debye- Waller factor is due to the scattering 
caused by heterophase fluctuations. 

15 Conclusion 

The systems are addressed, exhibiting phase transitions between thermodynamic phases with 
different symmetry, in which spontaneous symmetry breaking can be accompanied by local 
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spontaneous symmetry restoration, caused by the appearance of heterophase fluctuations. 
In such systems, above the phase transition point, where a disordered phase dominates, 
there appear the germs of the ordered phase. And below the transition point, where an 
ordered phase prevails, there arise the embryos of the disordered phase. Thus, around the 
phase transition point there is a region, where the phases with higher and lower symmetries 
coexist. The embryonic regions of a competing phase are distributed randomly in space. 
Their characteristic sizes are mesoscopic, so that the typical size r/ of such an embryo is 
much larger than the mean interparticle distance a, but much smaller than the characteristic 
size of the whole system L, that is, 

a < 77 < L . 

The typical size ry should be understood as an average size, since the shapes of the embryonic 
heterophase fluctuations are not necessarily regular, but can be rather ramified. Therefore, 
terming these germs mesoscopic could be better applied to the average number of particles 
Nf composing each of such embryos, as compared to the total number iV of particles in the 
system, so that 

K N f < N . 

This type of heterophase fluctuations is common for a number of substances. 

The general approach for treating these heterophase systems is developed, being based on 
the notion of weighted Hilbert spaces. The real-space distribution of the phases is described 
by means of manifold indicators. The averaging over random phase configurations reduces 
the problem to the consideration of an effective renormalized Hamiltonian composed of 
the phase-replica Hamiltonians representing the phases of different symmetry. Stability 
conditions define the geometric phase probabilities in a self-consistent way. 

The method is illustrated by applying it to heterophase ferroelectrics that are the typi- 
cal materials exhibiting the appearance of such heterophase fluctuations around their phase 
transition points between the paraelectric and ferroelectric phases. The influence of the het- 
erophase fluctuations is the strongest in the vicinity of the phase transition point. Numerical 
calculations show that the occurrence of such fluctuations leads to the noticeable decrease 
of the sound velocity and Debye-Walle factor at the transition point. The relative values of 
this decrease turn out to be universal, only weakly depending on the material parameters. 

In conclusion, it is worth mentioning that thermodynamic phases and phase transitions 
between them can be conveniently characterized not only by order parameters but also 
by order indices [HE]- Another important characteristic is the measure of entanglement in 
the considered physical system [97]. These three characteristics are interrelated with each 
other [981 - 1100] . The usual situation is when the increasing order is accompanied by the 
decreasing entanglement [IOOJ. Since a mesoscopic mixture is a system that is between 
an absolutely disordered and a completely ordered phases, its entanglement should be be- 
tween these two limiting cases. An additional entanglement arises in mesoscopic mixtures 
because a heterophase system consists of several spatially separated regions with different 
symmetry, and these mesoscopic regions are mutually entangled. This increases the system 
entanglement, as compared to the completely ordered case. The problem of studying the 
entanglement level of mesoscopic mixtures is extremely interesting. Since the mixture prop- 
erties can be regulated by means of thermodynamic parameters, this gives an additional 
possibility of governing the level of entanglement in such mesoscopic mixtures, which does 
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not exist in the case of pure systems. This novel possibility of regulating entanglement is 
really exciting, but it is a separate problem that goes out of the scope of the present paper. 
It is clear that the mesoscopic mixtures provide a nontrivial possibility for regulating the 
system entanglement, which can be used for quantum information processing. 
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